Macroscopic Quantum Coherence in Molecular Magnets 
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We study macroscopic quantum coherence in antiferromagnetic molecular magnets in the presence 
of magnetic fields. Such fields generate artificial tunnel barriers with externally tunable strength. 
We give detailed semi- classical predictions for the tunnel splitting in various regimes for low and high 
magnetic fields. We show that the tunneling dynamics of the Neel vector can be directly measured 
via the static magnetization and the specific heat. We also report on a new quantum phase arising 
from fiuctuations. The analytic results are complemented by numerical simulations. 
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Quantum spin dynamics in mesoscopic magnets has re- 
ceived much attention over the recent years, both from 
experiment and from theory ^ . A number of nanosized 
particles in the superparamagnetic regime have been 
identified as promising candidates for the observation 
of macroscopic quantum phenomena (MQP) such as the 
tunneling of the magnetization out of a metastable poten- 
tial minimum, or, more strikingly, macroscopic quantum 
coherence (MQC), where the magnetization (or the Neel 
vector) tunnels coherently between classically degenerate 
directions over many periods. On one hand, these phe- 
nomena are interesting from a fundamental point of view 
as they extend our understanding of the transition from 
quantum to classical behavior. On the other hand, the 
measurement of MQP quantities such as the tunnel split- 
ting provides independent information about microscopic 
parameters such as anisotropics and exchange constants. 

A prominent example of such MQC behavior that has 
attracted wide attention (but also scrutiny) is the anti- 
ferromagnetic ferritin . More recently, molecular mag- 
nets Q such as the ferric wheel or Mni2 have emerged 
as promising candidates for the experimental observation 
of MQP mainly for three reasons. First, molecular 
magnets have well-defined structures and magnetic prop- 
erties. Thus, precise values for the tunneling rates can be 
calculated. Second, molecular magnets can be produced 
as single crystals that contain a macroscopic number of 
identical magnetic subunits, which leads to a natural am- 
plification of the single-unit signal. Third, the typically 
high symmetry of these magnets reduces the number of 
independent parameters. 

In this letter we discuss novel tunneling scenarios in 
antiferromagnetic (AFM) molecular magnets. A key fea- 
ture of our discussion is to exploit the well-known fact ||] 
that an effective anisotropy can be generated in an AFM 
by applying a magnetic field. Thus it is possible to create 
tunnel barriers that are tunable by an external parame- 
ter. Evidently, such control parameters are highly desir- 
able as they open the door to systematic tests of MQC. 
We concentrate on ring- like structures such as Feg, Feio, 



Vs Q , where the spins interact with their nearest neigh- 
bors via exchange coupling. In particular, we show that 
the tunneling rates become field dependent and thus can 
be measured via the static magnetization and (less sur- 
prisingly) also via the Schottky anomaly of the specific 
heat. Our discussion is based on the non- linear sigma 
model (NLsM) that includes anisotropics and magnetic 
fields. The quantum dynamics of the Neel vector is then 
studied by instanton methods. Such methods are semi- 
classical in nature, i.e. valid for large spins and in the 
continuum limit. To cover the small end of the size scale, 
we performed ab initio numerical calculations; they agree 
well with the analytic results in the regime where a com- 
parison is possible. We find several distinct tunneling 
regimes, depending on the ratio of crystalline anisotropy 
to magnetic field. Motivated by recent measurements 
on single-crystal Fcio which indicate the presence of an 
anisotropy axis [Q, we give estimates of these MQC cor- 
rections in the magnetization and the specific heat, and 
we show that they are within experimental reach. 

We consider a ring-like molecular magnet, modeled as 
N spins s regularly spaced on a circle lying in the xy- 
plane, with N even. The Hamiltonian is (Sat+i = Si) 



(1) 



with AFM exchange coupling J > 0, and where Ui{Si) is 
the crystalline anisotropy at site z, h = 7B, with B being 
the magnetic field, 7 = g/is/Ti, and g is the electronic g- 
factor. For simplicity, we assume that both g and J are 
isotropic, and that the point symmetry of the molecule is 
that of a ring. Up to second-order in the spin variables, 
the most general form of the anisotropics is then Ui(Si) = 
kzSf^ + kr{Si ■ e^)^, where and kr are the axial and 
radial anisotropics, respectively, and is a unit vector 
at site i pointing radially outwards. We also assume that 
J 3> kz, kr, which is typically the case. 

We now derive an effective Lagrangian (NLsM) de- 
scribing the low-energy physics of (^ by extending stan- 
dard techniques |^ to include magnetic fields. We in- 
troduce spin coherent-states, and decompose the local 
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spin as Si = (— l)*sn + 1;, where the Neel vector n (with 
= 1) is taken as uniform for our small system, and 1^ is 
the fluctuation at site i (with Ij-n = 0). After integrating 
out the li , and keeping only the lowest-order terms in 
kz/J and k^/ J, we obtain the Euclidean Lagrangian 



[-{in A h - hf + (h-n)^] + NkzS^nl, (2) 



with a single effective axial anisotropy kz — kz — kr/2. 
Note that the magnetic field has two effects. First, it 
creates a hard axis anisotropy along its direction. This 
was already noticed in the context of spin chains Q , and 
is easy to interpret. The spins can gain Zeeman energy 
by canting towards the magnetic field, and this effect is 
maximal when the Neel vector is perpendicular to the 
field. What makes this anisotropy interesting for our 
purposes is that it is tunable by an external field. Second, 
a phase factor arises from the cross term 2ih • (n A ri). 
This is a crucial difference to the results of Ref . Q , which 
will have important consequences at high fields. 

Depending on the sign of kz and on the orientation of 
the field, various scenarios can be envisaged . For lack 
of space we restrict ourselves to the most interesting case 
where the field is in the ring plane, B = (5^,0,0), and 
perpendicular to a hard axis, i.e. kz > 0. The potential 
energy has then two minima at n = ±ej^. Tunneling 
of the Neel vector between these classically degenerate 
states results in a tunnel splitting of the ground state 
energy, which can be calculated using instanton methods 
pl| . For the semi-classical dynamics two kinds of hard 
axis anisotropies compete, the crystalline one, NkzS^, 
and that induced by the field, N{giJ,BBx)'^ /8J. Let us 
introduce their ratio A = 8s^Jkz/{giJ,BBx)^- For low 
fields (A » 1), the hard-axis is the z-axis, and the Neel 
vector, staying close to the xy-plane, tunnels via the x- 
axis. For high fields (A ^ 1) on the other hand, the 
hard-axis is the x-axis, and the Neel vector tunnels via 
the z-axis while staying in the j/z-plane. Without the 
phase term 2zh • (n A n) in (|^), the crossover would occur 
for A = 1, i.e. for a field B^ = s^/8Jkz/gfJ.B- As we will 
see, this extra term reduces the critical field. 

We first concentrate on the high- field regime (A <C 1). 
In this case, the Neel vector is conveniently parameter- 
ized as n = (cos 0, sin 6* sin 0, sin cos 0). We then find 
that the instanton solution minimizing the action belong- 
ing to (^) moves in the yz-plane with the frequency 
Whf = s\/8Jkz/h, and action 



Si,f/h = Nsy/2kz/J ± iTT 



Ng^BBx 
4J 



(3) 



where the upper and lower signs correspond to instan- 
tons and anti-instantons, respectively. The phase term 
in Eq. (^ arises because the spins cant towards the field, 
thereby acquiring an additional geometric phase factor. 

For the calculation of the fluctuation determinants, it 
is convenient to pass to dimensionless variables in (0). 



We write the action as S* = TiNgfiBBx/8J Jdr (L — h^), 
with 

hivif = gfiBBx, L = + 2i{n Ah)x + nl + Xnl, (4) 

where time is rescaled as f = luht. Expanding around the 
instanton solution up to second-order, one sees that the 
6- and (^-fluctuations are decoupled. For the ^-fluctuation 
determinant we find in leading order exp{±z-|-|-0(A'^/^)}. 
This is a new quantum phase distinguishing instantons 
from anti-instantons. Its occurrence is surprising, be- 
cause it is due to quantum fluctuations, and not due to 
a topological term in the action, in marked contrast to 
phases arising usually in spin problems Moreover, 
this new quantum phase does not depend on the spin 
s. The fluctuation determinant for (j) is standard [pT| , 
and by summing over all configurations |]l3|| , we finally 
find for the tunnel splitting of the ground state in the 
high-field regime. 
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Note that the tunnel splitting oscillates with the B-field 
as a result of interference between quantum spin phases 
(the new additional phase induces a shift from the usual 
cosine to a sine). Naively, one would have expected 
no field dependence, since both the tunneling barrier and 
the attempt frequency are constant. All these features 
are nicely confirmed by independent numerics (see below 
and Fig. 2). 

Next, we consider the low-field regime (A ^ 1), where 
we use the parametrization n = (sin cos 0, sin sine/), 
cos 9). We start by integrating out the (^-fluctuations 
to obtain an effective action. The term (n A n)^ in the 
Lagrangian (^ forces the instantons out of the xy-plane. 
However, for A ^ 1 the deviations are small. Thus, we 
write 9 — Tr/2 + 'd, and expand the Lagrangian to second- 
order in This gives L k, cfP' + cos^ (j) + 4ii?(/) cos (j) -\- "dG-d, 



where G = — 9| 4- A 



cos 



Integrating out is 



now straightforward. For A » 1 we find that G ~ X. 
Hence, we end up with the effective Lagrangian 



ieff - (1 
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(6) 



A simple quadrature shows that the instantons of this 
Lagrangian have an action SqB = 4(1 + 4/3A + 0(A^^)). 
Reinstating full units, we finally get for the action in the 
low-field regime (neglecting corrections of order A~^) 



2J 
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6 sV/c^ 



(7) 



Comparison with Eq. (^ shows that the crossover be- 
tween the low- and high-field regimes occurs for a mag- 
netic field B^ = asy/8Jk^/gfiB, with a = + 
(3 + VW)-^/^]/2 « 0.64, a sizeable reduction over the 
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result that would follow from neglecting the phase term 
2ih- (nAh) in (|). 

Next we determine the fluctuation-determinant. This 
raises one problem. While the action for the instanton 
solutions of is easy to calculate, the solutions them- 
selves cannot be obtained analytically. However, being 
only interested in the regime A S> 1, we can approximate 
this determinant by the one obtained for the fluctuations 
around the instantons of the Lagrangian L — (j)'^ + cos^ (p. 
This results in 
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The low-field splitting decreases (roughly) exponentially 
with the field. This is easily interpreted: The tunneling 
barrier increases quadratically with the field, whereas the 
attempt frequency increases linearly. 

We complete our derivation by discussing the range 
of validity of our results. A necessary condition to have 
tunneling (in the ground state) is that the barrier A[/ 
be much larger than half the attempt frequency uj p5| . 
Application of this criterion to the low- and high-field 
regimes yields two conditions, g^isBx 3> AJ/N, and 
NsyJkz/2J ^ 1. The effective Lagrangian (^) was de- 
rived under the assumption of (local) Neel order. Hence, 
the Zeeman energy must be smaller than the exchange 
energy, i.e. gfisBx ^ 4Js. Finally, our expression for 
high- field tunnel splitting is valid for A ^ 1, i.e. for 
g^sBx 2> sy^SJkzi whereas the low-field predictions 
hold if A 1, i.e. for g^sBx ^ sy/SJk^. We summarize 
m Fig. g the various regimes and critical fields we have 
obtained. 

We now turn to the question of how to observe the 
tunnel splitting. In contrast to previous cases jl^ such 
as ferritin ||] it is not possible to observe the switching 
of the Neel vector via an excess spin, since even if such 
an excess moment were present, it is easily seen |lo|] that 
it would always point along the magnetic field and not 
along the Neel vector. However, the dynamics of the Neel 
vector could be observed via resonances (occurring at A) 
in the NMR spectrum, which provides local spin informa- 
tion. An entirely different approach, which is only pos- 
sible because the tunnel splitting is B-field dependent, is 
to measure the static magnetization M — —l^g^B ^i) 
as a function of applied field. Indeed, we have seen that 
the two lowest energies are tunnel split by A and are sep- 
arated from the other levels by an energy huj ^ A, where 
A and lo are approximated by Aif, Ahf, and wif, Whf, re- 
spectively, depending on the field. At low temperatures, 
such that ksT <^ hw, the magnetization along the x-axis 
is then found to be 



, N 1\ A' / A 

Mx = [ ^gf^BBx ^- g^^B + - tanh ( ^ 



(9) 



Note that the first two terms in Eq. (p|) only give a linear 



dependence on Bx- The last term shows that deviations 
are proportional to A' = dA/dBx- 

The Lagrangian (|^) is not limited to the tunneling 
regime, but also covers the nearly- free limit of small 
kz, such that NsyJkz/^J <^ 1. This regime is most 
conveniently studied jl^] in terms of the correspond- 
ing Hamiltonian which is of rigid rotor type, -ffrot — 
-||^L^ -I- 7L • B -I- NkzS^nl, where n and the angu- 
lar momentum L satisfy standard commutation relations 
[Lj,nk] — ihejkini. For kz ~ 0, the ground-state is the 
state |/,/), with I = lgfiBBxN/4J\. Hence, the mag- 
netization Mx — IgfJ-B consists of steps of height g/is, 
occurring with a period gfis^Bx — iJ/N |18[| . This 
agrees with previous results obtained in the absence of 
anisotropies and tunneling [||. A small value of kz (be- 
fore tunneling sets in) leads to a rounding of the steps, 
as is easily seen from perturbation theory. For larger kz 
such that A ^ Tiu! the Neel vector no longer freely ro- 
tates but becomes strongly localized along the easy axis. 
As a consequence, the steps in the magnetization vanish, 
see (^). Conversely, notice that whatever the value of 
kz J sharp steps are recovered if the magnetic field is 
applied along the hard-axis. We have also confirmed this 
picture by direct numerical diagonalization of ffrot- 

The tunnel splitting is also reflected in other thermo- 
dynamic quantities. For example, the low temperature 
specific heat exhibits a characteristic Schottky anomaly 
|l9|, cv = kB{A/2kBTfsech^{A/2kBT), with a peak 
of height 0.64 /cs at a temperature T « O.GA/fcs. The 
location of this peak thus gives the tunnel splitting. 

The semi-classical analysis presented so far applies 
strictly speaking only to a sizable number of spins with 
s 3> 1. However, as is often the case with such methods 
the results are valid (at least qualitatively) even down 
to a few spins of small size. This expectation is indeed 
confirmed by direct numerical simulations which we have 
performed on Hamiltonian (|l|). Results for = 4 and 
s — 5/2, and for some typical values kz = J/10, and 
kr — Q, are presented in Fig. ||. We note that since most 
symmetries are broken in (|l|), larger system sizes become 
quickly inaccessible to numerics. The agreement with the 
semi-classical prediction is satisfactory in the high-field 
regime. Since for our test system the low- field regime 
becomes vanishingly small, i.e. 1 <C g^BBx/ J ^ •v/S, we 
cannot expect to find good quantitative agreement in this 
regime. Still, we can see from Fig. || that at the qualita- 
tive level the numerical and semiclassical approach show 
reasonable agreement. We have also calculated numeri- 
cally the matrix elements of the staggered magnetization 
and could confirm the tunneling picture. Thus, our the- 
oretical predictions give reasonably good results even for 
a very small cluster (and similarly for rings with larger 
N but smaller s |10|). Obviously, the accuracy of the 
semiclassical results will improve for larger systems. 

To support the experimental relevance of our results, 
we give some estimates for the ferric wheel, Feio, for 
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which N =10, s = 5/2, J/gnB = lOT §. While mag- 
netization measurements have been reported , no con- 
clusive comparison with our theory is possible presently 
since they have been performed on polycrystalline sam- 
ples with random orientation of the anisotropy axis, 
whereas the tunneling effect discussed here requires the 
B-field to have a fixed orientation with respect to such an 
axis. However, from the well-defined steps that have been 
observed and from recent single-crystal measurements 
one can infer that the magnitude olkz/ J is small, and of 
the order of 0.03 0. The low-field regime extends then 
from 4T to 7.8 T, with a tunnel splitting Ait//i decaying 
exponentially from (roughly) 6 • 10^° Hz to 4 • 10^ Hz. Cor- 
respondingly, the Schottky peak of the specific heat shifts 
from 1.6 to 0.12 K, while the tunneling corrections in the 
magnetization range from 60 % down to 16% of ^b- The 
high-field regime starts at 12 T, with the tunnel splitting 
having oscillations of magnitude Ahf//i ~ 6 • 10^ Hz and 
period 4T. The tunneling corrections in the magnetiza- 
tion reach at their peak 17% of jiB- The Schottky peak 
oscillates between zero and 0.2 K. The crossover temper- 
ature to the quantum regime, = fioj/AkB, is in the 
1-4 K range. Finally we remark that the same numbers 
apply to the high- field regime of an easy-axis system . 
Hence, all quantities appear to be well within experimen- 
tal reach. 
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FIG. 1. Summary of the various regimes we have obtained 
as functions of the applied magnetic field (a ~ 0.64). 
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FIG. 2. Results for the energy splitting Ai?i2 between the 
lowest two states (upper part) and the magnetization Mx at 
T = (lower part) for a system of A = 4 spins s = 5/2 with 
kz = J/10 and kr = 0. The symbols give the results of a 
direct diagonalization of the Hamiltonian (0), the continuous 
lines give the semi-classical predictions. The dashed vertical 
hues indicate the various critical fields (see Fig. m). 
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